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Abstract 

We study pilgrim dark energy model by taking IR cut-offs as par¬ 
ticle and event horizons as well as conformal age of the universe. We 
derive evolution equations for fractional energy density and equation 
of state parameters for pilgrim dark energy. The phantom cosmic evo¬ 
lution is established in these scenarios which is well supported by the 
cosmological parameters such as deceleration parameter, statefinder 
parameters and phase space of lo ^ and We conclude that the con¬ 
sistent value of parameter is < 0 in accordance with the current 
Planck and WMAP9 results. 
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1 Introduction 

Over the past decade, the substantial progress in astronomical observations 
indicate that our universe is presently going through the phase of acceler¬ 
ated expansion. Observations of type la supernovae (SNela) (Perlmutter et 
al. 1999; Riess et ah 2007), anisotropy measurement in current cosmic mi¬ 
crowave background (CMB) from WMAP (Spergel et al. 2004) and data of 
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large scale structure (LSS) from Salon Digital Sky Survey (SDSS) (Tegmark 
et al. 2004) strongly endorse this manifestation. The mechanism behind 
the expanding paradigm is usually assigned to exotic energy component with 
strong negative pressure entitled as dark energy (DE). This may determine 
the ultimate future of the universe but its cosmological origin and character¬ 
istics are still a complicated story. The most likely theoretical campaigner 
of DE is the cosmological constant A with equation of state (EoS) = —1 
(Weinberg 1989; Peebles and Ratra 2003). The model comprising of A and 
cold dark matter (CDM) dubbed as ACDM model which suffers from fine 
tuning and cosmic coincidence puzzles. However, the predictions of ACDM 
model appear to fit the current observational data (Samushia and Ratra 2008; 
Jassal et al. 2008). Numerous candidates of DE have been suggested in liter¬ 
ature such as quintessence > — 1), phantom {oj^ < —1) violating the null 
energy condition, quintom with evolving across —1, K-essence, tachyon, 
ghost condensate, holographic DE (HDE) and so forth (Arkani-Hamed et 
al. 2002; Armendariz-Picon et al. 2000; Caldwell 2002; Feng et al. 2005; 
Hsu 2004; Li 2004; Steinhardt et al. 1999). Introducing new ingredients 
of DE to the whole cosmic energy is one way to handle the issue of cosmic 
acceleration. Another approach is the modihcation of Einstein Lagrangian 
to get modihed theories such as f{R) (Sotiriou and Faraoni 2010), f{R,T) 
(Harko et al. 2011; Sharif and Zubair 2012a, 2012b, 2013a, 2013b, 2013c) 
and (Haghani et al. 2013; Sharif and Zubair 2013d, 2013e) 

gravities, T is the trace of the energy-momentum tensor. 

Cohen et al. (1999) set up a relation between ultraviolet (UV) and in¬ 
frared (IR) cut-offs due to the limit made by the formation of a BH. If 
is the quantum zero-point energy density associated with UV cut-off then 
entire energy in a sysytem of size L should not exceed BH mass of the same 
size so that ^ = I/VSttG is the reduced Planck mass. The 

largest IR cut-off saturates the inequality and one gets the HDE density 

3c2m2 

where 3c^ is a nnmerical constant. Several proposals have been suggested for 
IR cut-off including Hnbble, particle and event horizons as well as conformal 
age of the nniverse, Ricci scalar and Granda Oliveros cnt-off (Li 2004; Gao 
et al. 2009; Granda and Oliveros 2008, 2009; Wei and Cai 2008a). 

Phantom form of DE {p + p < 0) possesses a peculiar feature of big 
rip, the innumerous cosmic expansion within hnite time. In such scenario. 
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the energy density grows quickly and disrupts all the large structures and 
bounded objects. A question arises about the fate of BHs in the universe 
dominated by phantom DE. One can say that repulsive force would be strong 
enough to avoid the gravitational collapse and formation of BHs. Babichev 
et al. (2004) explored phantom energy accretion of BH and found that its 
mass decreases gradually. Some authors (Jamil and Qadir 2011; Sharif and 
Abbas 2011, 2012) discussed this issue for different BH solutions. Gao et al. 
(2008) showed that physical mass of BH may rather increase due to accretion 
of phantom energy implying the violation of cosmic censorship conjecture. 

Gonzalez and Guzman (2009) tested the accretion of phantom scalar field 
into BH with different initial conhgurations and found that this mechanism 
can reduce one half of the BH area. Sun (2009) studied dynamical equation 
of BH mass in terms of cosmological parameters and obtained that BH mass 
reduces to zero for the phantom dominated universe approaching to big rip. 
Recently, Wei (2012) proposed a new model of DE named as pilgrim DE 
(PDE) based on the idea that phantom DE is strong enough to avoid the 
formation of BH. He considered Hubble horizon as an IR cut-off and devel¬ 
oped constraints on PDE using the latest cosmic observations. Sharif and 
Jawad (2013) analyzed the interacting PDE models in terms of present day 
values of cosmographic parameters. 

This paper explores the cosmological evolution of PDE for three cut-offs 
namely particle horizon, event horizon and conformal age of the universe in 
FRW universe. We follow the work of Li (2004) in HDE to explore the fea¬ 
tures of non-interacting PDE for these cut-offs through fractional DE density 
EoS parameter statehnder diagnostic parameters and — analy¬ 
sis. The paper has the following format. In next section, we comprehensively 
present the evolutionary paradigm of PDE. We conclude our results in the 
last section. 


2 Pilgrim Dark Energy 

The pilgrim dark energy is dehned through the relation Wei (2012) 

( 1 ) 

The hrst Friedmann equation is given by 

= p, (2) 
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where p = pm + P-d comprises of matter as well as DE components and 
H is the Hubble parameter. The matter energy density is defined as pm = 
Pmo&~^^ {x = In a) from the matter energy conservation equation. By setting 
the fractional energy densities of matter and DE 

Dm = — , = —, Pcri = ^MlH\ 

Peri Peri 

Eq.(|2]) can be cast to the form 



(3) 


(4) 


If is known then one can determine the whole expansion history H{x). 
We discuss cosmological evolution for different cut-offs such as particle and 
event horizons as well as conformal age of the universe. 


2.1 Particle Horizon 


This horizon was initially used by Fischler and Susskind (1998) in holographic 
cosmology. Li (2004) discussed HDE by taking particle horizon as an IR cut¬ 
off and found that it does not imply realistic cosmology with EoS > —1/3. 
The particle horizon is defined as 


L 


Rp = a{t) 



a{t) 



(5) 


Combining the dehnition of PDE ([T]) and particle horizon (|5]), it follows that 



r da' 

Jo ^ 


1 

a V ) 


( 6 ) 


Equation ([3]) can be represented as 


1 _ o(l — D^) 

Ha HQ^/Vtmo 
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Substituting this relation in Eq.dH]), we have 



\/ a(l — 


na 


g(3//i-l)a: 


vai 



Differentiating it with respect to x = In a, it follows that 


% = D^(l - (3 - /i - ^(1 - D^)l/2-l/M^^^)l/Mg(3/2-3/M)x j ^ ^ 7 ) 


where C 


rj-2 —/j,q1 — m/2 
-^0 


1/m 


and prime indicates derivative with respect to 


X = In a. This result can explain the cosmic evolution according to PDE 
with particle horizon. 

One can exactly solve the above equation to represent the behavior of 
PDE (Li 2004). The corresponding EoS parameter can be set by using the 
energy conservation equation as 




Idhip^ 
3 dx 



1 / 

H J 


( 8 ) 


If we set /i = 2, Eqs.([7]) and ([8]) reproduce the corresponding results in HDE 
with particle horizon as an IR cut-off (Li 2004). For HDE, EoS is = 

— 1 -|- ^\/n^ indicating that when D,? —)■ 1 in the future, cu,? = —1 -f- ^ ^ 
which appears to be inconsistent with the accelerating phase. In case of PDE, 
we have dependence on parameter p resulting in extra degrees of freedom. 
As is involved in Eqs.d?]) and ([8]), so the expansion history is independent 
of signature of n. Solving Eq.(I7]) with the initial condition = 1 ~ ^mo 
and using in Eq.(|H]), the evolution of PDE is shown in Figures 1-4. Here, 
present day values of and H are defined from recent Planck results as 
Dmo = 0.315 and Ho = 67.3. The Planck and WMAP9 observations set the 
constraints for EoS of DE as ui^ = —1.131qjq and —1.71 < < —0.34 

respectively (Ade et al. 2013) 

For /i ^ 3, we have purely matted dominated phase of the universe as 
shown in Figure 1. This choice should be neglected in search of some ob- 
servationally consistent models. If /r < 0, the evolution of and are 
shown in Figure 2 which indicates that a;,? is always in phantom region and 
never intersects the phantom divide line = —1) in entire cosmic evolu¬ 
tion. This behavior is similar to the case of PDE with Hubble horizon in 
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Figure 1: Evolution of and oo^ for PDE with particle horizon. Here we 
set /i = 3 and n = 2. It clearly shows that such choice is no more realistic 
resulting matter dominated regime in the future evolution. 


which < —1 in whole cosmic history (Wei 2012). For the Hubble horizon, 
asymptotically goes to —1, he., it represents de Sitter phase in late times 
whereas in our case for PDE with particle horizon, it ends up with phantom 
phase. 

The deceleration parameter is dehned in terms of and D,? as 

aa 1 o \ 

Q = = -(1 + 

2 

which is a handy tool to explain the transition from decelerated phase to 
accelerating regime. We plot q versus z; and show the transition from decel¬ 
erated phase to accelerated era. The universe entered in accelerated era in 
recent past and it will hnish with q < —1 representing the phantom evolu¬ 
tion. The change of signature in q depends upon the values of p and the era 
of accelerated expansion begins earlier for large values of /i. Differentiating 
Eq. ([8]) with respect to x = In a, we get 

Caldwell and Linder (2005) discussed the quintessence feature of DE can¬ 
didate and analyzed its representation in plane. They established 

the limits of quintessence model in phase space of and and pointed 
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Figure 2: Evolution of and u:^ for PDE with particle horizon (p ^ —1). 


out two regions of this plane namely thawing > 0 with < 0) and 
freezing < 0 with uj^ < 0). It is remarked that cosmic expansion is 
accelerated in freezing region when compared with thawing region. This ap¬ 
proach has been applied in different settings by considering various forms of 
DE such as quintessence, phantom and quintom models (Chiba 2006; Gao 
et ah 2006; Scherrer 2006). In (Sharif and Zubair 2013f), we have also dis¬ 
cussed the phase space of a;,? and for new HDE which exhibits ACDM 
model (cu,? = — 1 and = 0) in future evolution. The evolution of in 
plane is shown in the right panel of Figure 3. This represents the 
freezing region for noninteracting PDE with particle horizon which favors the 
phantom evolution in this format of DE. 

In what follows, we examine PDE with particle horizon using the statehnder 
diagnostic. The pair {r, s} of statehnder diagnostic parameters is dehned as 
(Sahni et ah 2003) 


'a (r — 1) 

^ “ 3(g- 1/2)’ 

r is also named as jerk parameter. Statehnder diagnostic depends on the 
scale factor, its derivatives a and deceleration parameter q to diherentiate 
the cosmic expansion on geometric grounds. Cosmological models can be 
diherentiated on the basis of statehnder diagnostic as it shows distinct tra¬ 
jectories corresponding to specihc models. For ACDM model, the statehnder 
parameters are hxed as (r, s) = (1, 0) and in case of CDM regime these cor¬ 
respond to (r, s) = (1,1). In r — s plane, the trajectories for quintessence 
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Figure 3: Evolution of q and — plane for PDE with particle horizon for 
/i ^ — 1. A sign flip of q indicates transition to accelerated expansion and 
— io'g plane indicates the freezing region. The dots represent present day 
values of parameters. 


and phantom he in the range (s > 0, r < 1) whereas for chaplygin gas these 
correspond to (s < 0, r > 1). The statefinder diagnostic parameters can be 
represented in terms of oo^ and as 

3 —1 

r = 1 - [ui'^ - 3n;^(l + cn^)], s = -— [ui'^ - 3n;^»(l + . (11) 

Z oUJ^ 

Figure 4 shows the evolution trajectories of PDE with particle horizon in 
r — s and q — r planes. In the left panel, the evolution trajectories favor 
the chaplygin gas model with s < 0 and r > 1. We also plot the evolution 
trajectories of the deceleration parameter in g —r plane. Our results are con¬ 
sistent with (Wu and Yu 2005, 2006) where authors performed the statehnder 
diagnostic for the phantom and quintom DE model. This shows that the non¬ 
interacting PDE with particle horizon favors the phantom regime which is 
the basic idea of this candidate. Thus, for the realistic model of PDE, one 
needs to set /r < 0 and this choice is well supported by the results of Planck 
and WMAP9 observations (Ade et ah 2013; Bennet 2012). Here we take 
p ^ — 1 for which PDE implies that < —1 which is supported by other 
cosmographic parameters as shown in Figures 3 and 4. 
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Figure 4: Evolution of statefinder diagnostic parameters for PDE with par¬ 
ticle horizon (/i ^ —1). 


2.2 Event Horizon 

The IR cutoff L (event horizon) is dehned as 

r r. ^ ^ di ^ ^ da' 

Employing the dehnition of PDE ([1]) and event horizon (IT^ . we obtain the 
dynamical equation of fractional density of DE as 

D:, = D^(l - (^3 - /r + ^(1 - j _ ^^3^ 

The time derivative of PDE with event horizon as an IR cut-off is 


P-& — L — HL — 1. 


Using the energy conservation equation of DE, we obtain 


oj^ — —1 -|- 


P 




H \n^Mp~^ 


-i( 


(14) 


One can reproduce the corresponding results in HDE with event horizon 
for /i = 2. In case of HDE with event horizon, the EoS parameter is uj^ = 


9 















(a) 



z 


(b) 


-1.0 






-1.5 




/ 

1 

-2.0 

■ — 

- - ,,=-1 



— r - T - 

-2.5 





- 

-3.0 


--- 11=3 

-3.5 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 


Z 


Figure 5: Evolution of and versus z for PDE with event horizon for 
n = 0.5 and different values of fx. The dots indicate the present day values. 

— I which can result in three signihcant eras of cosmic expansion. 

If the universe is dominated by DE components, he., —> 1 in the future 
then for n > 1, oj^ is always greater than —1 which depicts quintessence era 
so that the universe escapes from getting in de Sitter and big rip phases. For 
n = 1, the universe enters the de Sitter era in future evolution and n < 1 
represents phantom phase where the universe behaves as quintom model of 
DE as intersects the cosmological constant line. The value of n plays a 
vital role in deciding the evolutionary features of HDE and ultimate fate of 
the universe. 

For PDE, the role of p is more crucial as compared to that of n. In this 
setting, we are mainly concerned with the choice /i < 0 but for PDE with 
event horizon, one can also set /i = 3. The evolution trajectories of EoS 
and fractional energy density of PDE are shown in Figure 5. For /i ^ — 1, 

approaches to 1 as z ^ —1 showing that DE dominates in later times 
of the universe. The EoS parameter uj^ is in the phantom regime (Figure 
5(b)) and the present day values of uj^ are consistent with the Planck results 
showing oj^ = —l.lSlgjo (Bennet 2012). In case of /i > 2, we show the 
evolution trajectories for /i = 3 which does favor the phantom feature of DE. 
For fi = 3, the curves in Figure 5 show somewhat distinct behavior where 

< —3.5 as z ^ —1. We neglect the values of /i > 3 because these values 
do not imply realistic results. The evolution of q versus 2 ; is represented 
in Figure 6(a) which shows that the universe entered into the accelerating 
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Figure 6: Plot 6(a) shows the variation of q versus z and 6(b) represents 
evolution trajectories of u:^ — for non-interacting PDF with event horizon 
corresponding to n = 0.5 and different values of fi. The dots indicate the 
present day values of q and 


phase in the recent past and it would switch over to q < —1 indicating the 
phantom paradigm. 

Taking derivative of Eq. ffTTj) with respect to x, we get 

+ (15) 

The phase plane of and for non-interacting PDF with event horizon 
is shown in Figure 6(b). The plane — uj'^ represents the thawing region 
in the evolution of PDE with event horizon. In the limit of future evolution 
z —)• — 1, —)■ 0 and > —1. Figure 7 shows statehnder analysis in 

q — r plane for the choice p ^ —1. The evolution trajectories in g — r plane 
for PDE with event horizon start from (g < 0.5, r < 1) and end up with 
(g < —1, r > 1). 

2.3 Conformal Age of the Universe 

The two time scales, age of the universe and conformal time have been sug¬ 
gested in literature corresponding to agegraphic DE (ADE) (Cai 2007) and 
new agegraphic DE (NADE) (Wei and Cai 2008a). These models can derive 
the cosmic expansion consistent with the recent observational data (Wei and 
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Figure 7: The statefinder analysis of PDF with event horizon in g — r and 
r — s planes. 


Cai 2008b) which can resolve the causality problem. However, it is pointed 
out that ADE model is classically unstable and NADE is no better than 
HDE in explaining the DE dominated universe. The NADE has been stud¬ 
ied to address various cosmological issues in Einstein and modihed gravities 
(Karami 2010; Jamil and Saridakis 2010; Sheykhi 2010). The conformal age 
of the universe is dehned as 

Using Eqs.(II]) and flTOl) . the rate of change of fractional DE density is 


The time derivative of is obtained as 

The corresponding EoS parameter is 


UJ'ff — —1 + 


/i 


( 




3aH 


il/i 


( 18 ) 


To demonstrate the evolution trajectories for NADE version of PDE, we 
include some facts about non-interacting NADE model. Wei and Cai (2008b) 
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showed that coincidence problem can be resolved for NADE if one chooses 
the value of parameter n nearly unity. The NADE is constrained from the 
observational data of SNela, CMB and LSS which implies the best £t value 
ofn = 2.76«:l“ (with la uncertainty). They found that EoS parameter for 
NADE approaches to —1 in later times regardless of the value of n. Zhang et 
al. (2013) showed WMAP 7-years observations set appropriate measure of n 
as n = 2.673 ]'^q;q771q;}5[1q; 222- In previous study (Sharif and Zubair 2013f), we 
have reconstructed f{R) gravity corresponding to NADE and set the values 
of n as n = 2.3, 2.8, 3.3, 3.8. It has been shown that these parametric values 
support the ACDM model in future evolution. 

For PDE with conformal time scale, we are mainly concerned to explore 
the behavior of parameter /i in determining the evolution for conformal time. 
Initially, we set /r < 0 and found that /r = — 1 supports the cosmological 
constant regime in future evolution as shown in Figures 8 and 9. Figure 
8(a) shows that —)■ 1 in the future evolution so that our universe is 

dominated by DE. is always less than —1 in the whole cosmic history and 
it will asymptotically approach to —1 in the future evolution of the universe. 
The present behavior of EoS parameter favors the phantom DE consistent 
with Planck and WMAP9 results. Therefore, evolution of the universe will 
end up with ACDM model rather than big rip. Such behavior is identical 
to that suggested by Wei and Cai (2008b) for PDE with Hubble horizon. 
In Figure 9(a), the evolution of deceleration parameter shows the bouncing 
behavior of the universe which entered in the expansion phase in recent past 
and it would end up in a de Sitter phase. 

If we differentiate Eq. ffTSj) with respect to x, we have 

= + ( 19 ) 

Figure 9(b) shows the evolution trajectories of in — plane for PDE 
with cosmological time scale. It is obvious that evolution in phase space of 
and results in ACDM model = —1, = 0) as 2 ; —>• —1 (or x —)■ cxd). 

We analyze that if one sets —55 < fi < —1, it would result in a;,? —)■ 0 {i.e., 
the matter dominated universe) in future evolution. If one chooses fi ^ —55, 
we can have w,? > — 1 in later times of cosmic evolution but in such case the 
present day value is only consistent with WMAP9 observational results as 
it represents the quintessence model of DE (see Figure 10). For PDE with 
cosmological time scale, we also present the evolution trajectories for fi ^ 3 
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Figure 8: Evolution trajectories of and for PDE with conformal time. 
We set /i = — 1 and n = 2. 
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Figure 9: Evolution trajectories of q and uj^ — phase plane for /i = —1 
and n = 2. Star indicates the ACDM model with = —1 and u:'^ = 0. 
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Figure 10: Evolution trajectories of uj^ and q versus redshift for PDE with 
conformal time corresponding to fi = —55 and n = 2. 


as shown in Figure 11. It is found that for /x = 3, we have ACDM regime 
whereas for /x ^ 4, it results in quintessence era of the universe. 


3 Conclusions 

In this paper, we study the phantom evolution of PDE with three cut-offs 
namely particle horizon, event horizon and conformal age of the universe in 
FRW spacetime. We explore these IR cut-offs to establish the consistent 
range in PDE for parameter fi and also the phantom regime. Following 
(Li 2004), the evolution equation of fractional energy density of DE and 
dynamical relation of are formulated in these settings, li fi = 2, one 
can determine the corresponding results in HDE with particle horizon, event 
horizon and conformal age of the universe. We set the present day values of 
parameters according to recent Planck observations and present the evolution 
inconsistent with this data set. 

Firstly, we have analyzed the non-interacting PDE with particle horizon 
and shown the evolutionary paradigm of and If fi ^ 3, then we have 
purely matter dominated phase of the universe since —)■ 0 and > 1 for 
2 : —)■ — 1 as shown in Figure 1. This choice is neglected in search of some 
consistent models. For fi < 0, the plots of and are shown in Figure 
2 which represents that is in the phantom regime and never intersects 
the phantom divide line in the whole cosmic history. This is identical to 
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Figure 11: Evolution trajectories of and q versus 2 ; and in 

and q — r planes for PDE with conformal time. Plot (a) shows the evolution 
trajectories of which results in DE dominated era in future evolution. 
In plot (b) it is evident that for /r = 3, we have de Sitter phase whereas 
/i ^ 4 represents the quintessence regime. This behavior is well supported 
by the evolution of q in plot (c). The evolution trajectories of are shown 
in plot (d) which indicate the ACDM = (—1,0) model for /i = 3. 

Statefinder diagnosis can also be seen in plot (e). 
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that of PDE with Hubble horizon (Wei 2012). The evolution of shows the 
phantom regime inconsistent with the current observational results of Planck 
and WMAP9 data sets. We also plot deceleration parameter and the phase 
space of and in Figure 3. lies in the freezing region favoring the 
phantom evolution in this format of DE. The statehnder diagnostic (Figure 
4) show that evolution trajectories are consistent with phantom regime (Wu 
and Yu 2005, 2006). Thus, we conclude that for realistic model of PDE with 
particle horizon one needs to set p < 0 and this choice is well supported by 
the results of Planck and WMAP9 observations. 

Secondly, we have explored PDE scenario in the light of event horizon. 
We are mainly concerned with the choice of p < 0 but for PDE with event 
horizon one can also set p = 3. For p ^ —1 and p = 3, the evolution of 

and is shown in Figure 5. It is found that DE dominates in future 
evolution as —)■ 1 for approaching to —1. The EoS parameter is in the 
quintessence regime in recent past which bisects the phantom divide line and 
ends up in phantom era. For p > 2, the acceptable results are found only 
for fi = 3 whereas parameter p > 3 does not show realistic results. Figure 
6(a) favors the phantom DE showing sign flip of q in recent past leading to 
q < —1. We also show the evolution trajectories in plane (Figure 

6(b)). It represents the thawing region and in later times of the universe 

—)■ 0 with > —1. These results are also presented in q — r plane (Figure 
7). It is found that for PDE with event horizon the acceptable range of p is 
p ^ —1 and /i = 3. 

Thirdly, we have used conformal time scale as IR cut-off for PDE. Initially, 
we set p = —1 and plot the fractional density and EoS parameter of DE 
(Figure 8). Accordingly, D,? —)■ 1 shows the dominance of DE and < —1 
in entire cosmic evolution which approaches to —1 in ultimate fate of the 
universe. Consequently, cosmic evolution ends up with cosmological constant 
regime avoiding the big rip singularity which is identical to Wei (2012) for 
Hubble horizon. The deceleration parameter q shows the bouncing behavior 
of the universe where the universe entered in accelerated expansion era in 
recent past and concludes in de Sitter phase. These results are also favored 
by the phase space of and as shown in Figure 9(b). It is found that 
the range —55 ^ /i < —1 results in matter dominated cosmic evolution. If 
one sets p ^ —55, can show phantom evolution for later times but in such 
case the present day value of is not observationally consistent. Hence, for 
conformal age of the universe in PDE scenario, the only acceptable value of 
jj, is —1. 
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In (Sharif and Jawad 2013), anthers discussed the PDE for non-interacting 
case by dehning EoS and other cosmographic parameters in terms of present 
day values of Qm and if. This study is conhned to present scenario and 
does not show the behavior in entire cosmic evolution. Comparatively, the 
dynamical equation of with the initial condition = 1 ~ ^mo implies 
that the entire cosmic evolution of fi,? as well as cj^ can be established. We 
remark that PDE with particle horizon leads to phantom evolution if p < 0 
and for event horizon one can set both ju < 0 and /i = 3. In case of conformal 
age of the universe, the result is quite signiheant where only consistent value 
is /i = — 1 showing identical behavior to that for the Hubble horizon (Wei 
2012 ). 
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